Decoherence in a system of many two— level atoms 
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I show that the decoherence in a system of A'' degenerate two-level atoms interacting with a 
bosonic heat bath is for any number of atoms A*" governed by a generalized Hamming distance (called 
"decoherence metric" ) between the superposed quantum states, with a time-dependent metric tensor 
that is specific for the heat bath. The decoherence metric allows for the complete characterization of 
the decoherence of all possible superpositions of many-particle states, and can be applied to minimize 
the over-all decoherence in a quantum memory. For qubits which are far apart, the decoherence is 
given by a function describing single-qubit decoherence times the standard Hamming distance. I 
apply the theory to cold atoms in an optical lattice interacting with black body radiation. 



Decoherence is considered one of the main ingredients 
in the transition from quantum mechanics to classical 
mechanics [ll. According to experimental evidence (see 
for example 0, 01 ) "decoherence" can be well described 
within the standard framework of quantum mechanics 
01 . The loss of coherence of quantum mechanical super- 
positions becomes increasingly rapid when the "distance" 
between the components of the superpositions gets large 
0, 0, such that quantum mechanical behavior like in- 
terference of wave functions is never observed for macro- 
scopic objects. With few notable exceptions H, 
most work on decoherence has focused so far on systems 
with only one or a few effective degrees of freedom cou- 
pled to a heat bath. Even when large numbers of particles 
come into play, like in superconducting qubits based on 
Josephson junctions ^3 or superradiance 0,^3' '^^^y 
a single effective degree of freedom is normally retained. 
With the rise of quantum information science the funda- 
mental question whether for large particle (qubit) num- 
bers there rnight be new self-limitations of quantum me- 
chanics 0,0, 01, based e.g. on the amount of entan- 
glement in the quantum states 0, 0| , has also become 
of huge practical importance. Indeed, understanding and 
possibly suppressing decoherence of quantum superposi- 
tions in the exponentially large Hilbert space of possibly 
thousands of qubits is considered the most difficult prob- 
lem on the way to a large scale quantum computer. 

In the following I show that a generalized Hamming 
distance with a time-dependent metric tensor given by 
the heat bath and the couplings to the heat bath com- 
pletely determines the decoherence process of a an arbi- 
trarily large number of qubits. The standard Hamming 
distance, which counts the number of bits by which two 
code words differ, plays a central role in both classical 
and quantum error correction |l8l Hij , and determines in 
particular what code words remain distinguishable in a 
given error model. It is therefore very satisfying that the 
heat bath itself determines the metric relevant for the de- 
coherence process, as decoherence reflects to what extent 
the heat bath distinguishes the superposed states [HI. 

Starting point of the analysis are N two-level atoms 
at arbitrary but fixed positions {i — Q, . . . , N — 1) 



interacting with a common bosonic heat bath. All atoms 
are assumed identical with level spacing MIq and energy 
eigenstates | - 1) and with cr^] ± 1) = ±| ± 1). The 
total Hamiltonian reads 



(1) 



k i=0 

where axi and cr^i are Pauli matrices for atom i. While 
the results obtained are valid for a general bosonic heat 
bath, we will consider in the following as concrete exam- 
ple thermal black body (BBR) radiation dipole-coupled 
to the atoms |2^. The index k then stands for wave 
vector k and polarization direction A (kj — 2'inij/L 
with integer Uj, j = x,y,z for periodic boundary condi- 
tions) of the electro magnetic waves; (a^,) are the cre- 
ation (annihilation) operators for mode k with frequency 
LOk — c|k|, polarization vector efc, and electric field am- 
plitude Ek — \J fiajk/i'^SoV), where £q, c, and V are the 
dielectric constant of the vacuum, speed of light, and the 
quantization volume, respectively. The coupling constant 



of atom i to mode k is denoted by g\. 
where iS^^ stands for a unit vector in the direction of the 
dipole moment of atom i, (— l|d|l) — edu'-^^ with elec- 
tron charge e and dipole length d. Hamiltonians of the 
form 1^ have been studied before in many situations, in 
particular in the context of dipole-dipole interactions in 
quantum optics and various forms of mode selection (see 
[21| and references therein). We will restrict ourselves to 
atoms with degenerate energy levels, fig = 0. This should 
set a lower limit for decoherence in quantum memories 
based on qubits with an allowed dipole transition. From 
a theoretical perspective degenerate levels are extremely 
attractive, as the model can be solved exactly for any N. 
The results remain valid for small but finite ^Iq, as long 
as the time t <C l/(Nhflo) 0. The results are also eas- 
ily extended to any model with commuting system and 
interaction Hamiltonian, [i?sys, ^^int] = 0, such as pure 
dephasing models with i?sys CTzi and iJint c>c azi- 
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It is interesting to consider the decoherence process of 
n selected atoms (i = 0, . . . , n — 1) out of the N atoms. 
This offers the additional information how decoherence 
scales as function of the size of a sub-cluster in the quan- 
tum memory. The common heat bath can lead to entan- 
glement between the atoms li^,!!^!, such that unobserved 
atoms can become an important source of "indirect" de- 
coherence. Besides its conceptual importance, this can 
be relevant practically, if additional atoms are trapped 
in an optical lattice, or when information is stored in an 
atomic gas in a large number of atoms • We therefore 
first trace out the bosonic modes and secondly the ad- 
ditional atoms n . . . N — 1. The reduced density matrix 
p of the remaining atoms will be expressed in the eigen- 
basis of the Uxi, the natural basis (also called pointer 
basis) for studying the decoherence process j^. The ma- 
trix elements of p are pgg'(t), where s = (sq, . . . , Sji-i) 
is a subset of the code word s = (sq, si, . . . sn-i) (with 
Ca;i|sz)a; = Si\si)x, Si = ±1, and similarly for s', i.e. s and 
s' are binary representations of the matrix indices). We 
assume that the heat bath is initially in thermal equilib- 
rium at temperature T, and all atoms are prepared in a 
pure state. 

Using the method of shifted harmonic oscillators [l^ 
one shows that the time evolution of p is given by 



Pss'it) = Pss'(O) X 



(2) 



f n—1 
n-lj-l 

+i^^{s^Sj - s-s'j)ip^j{t,R^ - Rj) ^ 

j=l i=0 ■' 
N~-l /n-l \ 

n ~ s'i}tpii{t,Ri - R;) . 



l—n \i—0 J 

The functions and ipij are given by 

/y(t,R) = ^^i^cos(k-R)(l-coswfci)coth^^ 

(i) (i) 

(/7y(i,R) = 2V^^^^cos(k-R)(a;fct-sina;fet). (3) 

fe ""fc 

They both vanish at t = 0. 

We are interested in the initial stages of the decay 
of the off-diagonal matrix elements (the "coherences"), 
when jij and ipij are both much smaller than one, as 
later basically no coherences are left anyway. I define the 
"decoherences" dgg'{t) through 



\Ps,'{t)\ 



1 -dss'(t). 



(4) 



|Pss'(0)| 

Expanding eq.|(21) to lowest order in and ipij we find 



where denotes the transpose, and ]V[(t) is a real, sym- 
metric, and non-negative matrix with elements (i,j — 
0,...,n-l) 

My = 4/„ (i,R, -R,) + 2$y(i,R„Rj)(6) 
$y (t, Ri, Rj) = ipait, Rj - Ki)ipjiit, Rj - R,) . 

l—n 

The properties of M allow us to consider M(t) as a time- 
dependent metric tensor, and to define a metric in the 
vector space containing all code words. The heat- 
bath itself therefore induces a natural distance defined 



as 



|s-s-|U/(t) = -V(§-§')M(t)(i-§') 



(7) 



We will call ||s— s'||M(t) with M.{t) given by eq.® "de- 
coherence metric" and M "decoherence metric tensor" 
(DMT). Eq.© shows that the decoherence process of a 
superposition of code words s and s' is governed by the 
decoherence metric, 



l§ §'llM(t) 



(8) 



dss'it) ~ -(s-s')M(t)(s-s')^, 



(5) 



with a DMT M{t) defined by the heat bath itself and the 
couplings to the heat bath. Eq. ((SJ constitutes the central 
result of this letter. 

The distance ||s — s'||M(t) generalizes the well-known 
Hamming distance Z?^ (§,§'), which is defined as the 
number of digits in which s and s' differ. The Ham- 
ming distance D^{s,s') is recovered from eq.© for the 
trivial metric tensor M = I, where I is the identity ma- 
trix in n dimensions, as D^{s,s') = ||s — s'Wj. As an 
immediate consequence we find that the entanglement of 
a quantum state does not (or at least not fully) deter- 
mine the decoherence between its components. Indeed, 
all states related by parallel qubit flips in the components 
(such as e.g. (|000) + |111))/V3 and (|100) + |011))/V3 
for TV = 3) have the same entanglement and the same 
Hamming distance, but the off-diagonal elements of M 
lead to different time dependent decoherences for these 
states (see also FIG|21l. 

The decoherence metric is in the strict sense a pseudo- 
metric. It is symmetric, non-negative, and obeys the 
triangle inequality, but there can be code words s and 
s' with s 7^ s' such that ||s — s'||M(t) = 0. Such code 
words constitute a decoherence free subspace (DFS), 
which plays an irnportant role in quantum information 
processing |2a, 1221, |23, |22| • For the ease of language, I will 
continue to use the word "metric" instead of "pseudomet- 
ric" , however. The desire to keep the triangle inequality 
motivates the choice of the square root in eq. ^ , whereas 
the Hamming distance is defined without it. The latter 
still obeys the triangle inequality for binary code words, 
but the triangle inequality can be broken if the metric 
tensor has non-diagonal matrix elements and the met- 
ric is defined without the square root. The decoherence 
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metric has the form of a Mahalanobis distance, frequently 
used in statistical analysis t3_0| . 

The first part of the DMT, Af,j(t,Ri - Rj) in eq.®, 
depends on the n selected atoms only. It describes the 
direct decoherence of the atoms due to their exposure to 
the heat bath. The second part, 2$^ (t, R^, R^), however, 
depends on all atoms, even the non-selected ones, such 
that this contribution is not translationally invariant. $ij 
describes the "indirect decoherence" which arises due to 
the effective interaction and subsequent correlations or 
even entan glem ent that the heat bath induces between 
the atoms |22l Esj . Remarkably, the total decoherence 
dtotit), as defined by Es^^s' Ip^s' (01/ Es^^^s' Ifss' (0)1 = 
1 — dtot{t), involves only the trace of the DMT, dtot{t) — 
2_2-,.+i trM(t). In situations where is absent (e.g. for 
n = TV), the total decoherence scales for large n propor- 
tional to n, as fa is independent of n; otherwise, the 
scaling of with n and N has to be taken into account. 



The DMT can be evaluated explicitly for the present 
physical model. We will limit ourselves to direct deco- 
herence of n = two level atoms interacting with BBR, 
such that = 0, delegating indirect decoherence to an- 
other publication. We express all lengths in terms of the 
dipole length d and times in units of d/c, t ^ ct/d, and 
define — [Ri — Rj\/d. To simplify the presentation 
we will assume that all dipoles are oriented in the same 

(i) 

direction. This means that g/. is independent of i, and 
fij{t,'Ri — Rj) depends on i and j in the frame aligned 
with Ri — Rj only through and the angle i? between 

the Wj, and R^ — Rj. The azimuthal angle is irrelevant, 
and we can write /ij (i, Ri — Rj) — f{t,rij,'d). In the con- 
tinuum limit V ^ oo, the sums in eq.© become integrals 
over k which can be done analytically for T = 0. The re- 
sult reads f{t, r, ■(?) = (sit, r) sin^ d + c(t, r) cos^ i?) 
with 



1 - (t/r)^ 
sin nr 



cos(Kr) cos(Kt) H — sin(Kr) sm{Kt) — ( - 
r \ r 



c{t,r) 



+2 

2- I CifKlt - r\) - Ci(K|t -f r\) - In 
r 





t - r 








t + r 



t - r 



t + r 



- Gi(K\t ~ r\) + C\{K\t + r\) 



H sin(Kr)(cos(Kt) — 1) , 



(9) 
(10) 



where a = e^/(47reoSc) ~ 1/137.06 is the fine-structure 
constant, and k = kj^axd a UV cut-off of the integral 
over k. A cut-off is only needed for the diagonal part 
fa which corresponds to Vij = 0, but in order to assure 
the positivity of the decoherence metric, the same cut-off 
should be used for both diagonal and off-diagonal parts. 
The off-diagonal terms fij arise due to interference effects 
between two qubits. They decay as functions of r on a 
length scale of the order 1/k (see FIG^). There can also 
be regions in the t, r plane where fij becomes negative 
(not withstanding the fact that the decoherence metric is 
always non-negative) . One such zone is a ridge along the 
light-cone t = r. Another one exists for fixed r almost 
independently of t, for distances of the order 1/k. 

The function /(t, r, -d) is continuous in the limit r — )■ 0, 
and leads to diagonal matrix elements fa independent of 
•d and i, fii{t,0) = f{t,0,'d) = fii{t), which describe 
single-qubit decoherence (i.e. n = 1). Explicitly, 

2 f 1 — cos(Kt) — sin(Kt) \ 
^^^^')=3^"U^ ^ )- ^''^ 

The expressions H9I1()I11|I are exact in the limit T 
0. Corrections due to finite temperature are of order 
ksT/ (huji-aax) with Wmax = cfcmax, and will be neglected 



in the following. They can be easily taken into account 
by evaluating the integrals obtained from eqs.© numeri- 
cally. The divergence of fa for k oo makes the cut-off 
a relevant physical quantity. A cut-off always arises, at 
the latest at k ~ 1 when the dipole approximation breaks 
down. Much smaller values of k are expected for realistic 
atoms when transitions to excited levels come into play. 

It is clear from eq.© that in the limit of large qubit 
distances (r ^ 1/k) the decoherence between two code 
words is just given by the Hamming distance times the 
single-qubit decoherence, dss'(i) = 'ifoo{t)D^ {s,s'), as 
the qubits are then independent and each qubit con- 
tributes to decoherence if and only if it is in a super- 
position of I — 1) and |1). This is shown in FIG[21for the 
example of 3 x 3 qubits arranged on a 2D square optical 
lattice with lattice constant a. The larger a the smaller 
the spread of the decoherences around the value given by 
the Hamming distance and single qubit decoherence. 

The minimum and maximum decoherence, and the 
width of the decoherence distribution, show a substantial 
dependence on the lattice spacing (see FIG 13). Notably, 
for small lattice spacing (r <C 1/k) the minimum deco- 
herence tends to zero. A DFS arises, as the atoms then 
couple in a symmetric way to the heat bath (the latter 
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FIG. 1: (Color online) The function f{t,r,j,-d), for n = 0.01 
and i9 = 7r/2atr = 0. The limit r — > represents the 
single-qubit decoherence, fii{t). 
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FIG. 2: (Color online) Decoherences ds,s'(i) of 9 qubits in 
a 3 X 3 square optical lattice a,t t — 200.0 as function of 
the Hamming distances D^(s,s'). Small red squares are for 
a — 580d, large blue circles for a = lOOOd. Same parameters 
as in FIGIU 




a/d 



FIG. 3: (Color online) Statistical analysis of all 2^'' indepen- 
dent decoherences of 3 x 4 qubits in a 2D square lattice for 
a = lOOd . . . llOOd (parameters as in FIG|5J. The average de- 
coherence, given by the trace of the DMT, is independent of 
the lattice spacings (central green curve, error bars represent 
± one standard deviation). The minimum decoherence (bot- 
tom curve, blue) indicates a DFS at small spacings; the top 
curve (red) is the maximum decoherence. The inset shows a 
histogram of all decoherences for a = 580d. 



does not contain sufficiently short waves for distinguish- 
ing the positions of the atoms). The decoherence metric 
allows to generalize the DFS concept in the sense that 
even for non-symmetric coupling one may optimize now 
the performance of the quantum computer by encoding 
in the states with smallest decoherence, which can be 
explicitly found using the decoherence metric. In |3l| | a 
method was introduced that allows to quantum compute 
within a DFS, based on a strong coupling to the environ- 
ment and corresponding broadening of non-DFS states. 
It is expected that this method can be generalized to the 
current approach. 

As a summary, I have derived a "decoherence metric" 
that induces a natural distance between binary quan- 
tum code words. The decoherence metric generalizes the 
well-known Hamming distance and determines directly 
the time-dependent decoherence in an arbitrarily large 
system of qubits with degenerate energy levels coupled 
to a heat bath of harmonic osciallators. 
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